SPARTA (26 Feb 2021)
################################################################################
# beam of particles striking the surface at an inclined angle
# free molecular flow (no collisions)
#
# Note:
#  - The "comm/sort” option to the “global” command is used to match MPI runs.
# The "comm/sort" option should not be used for production runs.
################################################################################

seed	    	    123456
dimension   	    3
global              gridcut 0.0 comm/sort yes

boundary	    	oo oo so


create_box          -11 11 -11 11 0 10
Created orthogonal box = (-11 -11 0) to (11 11 10)
create_grid 	    2 2 2
Created 8 child grid cells
  CPU time = 0.00138879 secs
  create/ghost percent = 92.4635 7.53648
balance_grid        rcb cell
Balance grid migrated 0 cells
  CPU time = 0.000130653 secs
  reassign/sort/migrate/ghost percent = 53.1022 1.64234 31.2044 14.0511

global		    	nrho 1e10 fnum 1e6

species		    	air.species O CO CO2 O2 C
mixture		    	air O O2 vstream 0 1000 -1000

mixture             air O   frac 1.0
mixture             air CO  frac 0.0
mixture             air CO2 frac 0.0
mixture             air C   frac 0.0
mixture 			air O2 frac 0.0


surf_collide        1 cll 300.0 0.5 0.5 0.5 0.5

bound_modify 		zlo collide 1

##################################### SURF REACT ADSORB ######################################
##################################### FACE/BOUNDARY OPTION ###################################

#surf_react        	adsorb_test_ps1 adsorb ps sample-PS_1.surf nsync 1 face 1000 6.022e18 O CO
#bound_modify        zlo react adsorb_test_ps1


surf_react        	adsorb_test_ps2 adsorb ps sample-PS_2.surf nsync 1 face 1000 6.022e18 O CO
bound_modify        zlo react adsorb_test_ps2

########################## BEAM ############################################################
# Beam at multiple points so that different processors handle the surface collisions

region              circle1 cylinder z  0 -10 1 -INF INF

fix                 in1 emit/face/file air zhi data.beam beam_area_1 nevery 100 region circle1

################################################################################################

#dump                2 image all 10 image.*.ppm type type pdiam 0.2 surf proc 0.01 size 512 512 zoom 1.75 gline no 0.005
#dump_modify	     2 pad 4

timestep            0.0001

stats		    	10
stats_style	    	step cpu np nattempt ncoll nscoll nscheck
run 		    	1000
Memory usage per proc in Mbytes:
  particles (ave,min,max) = 0 0 0
  grid      (ave,min,max) = 1.51379 1.51379 1.51379
  surf      (ave,min,max) = 0 0 0
  total     (ave,min,max) = 1.51379 1.51379 1.51379
Step CPU Np Natt Ncoll Nscoll Nscheck 
       0            0        0        0        0        0        0 
      10 4.8875809e-05        0        0        0        0        0 
      20 0.00012159348        0        0        0        0        0 
      30 0.00020384789        0        0        0        0        0 
      40 0.00027489662        0        0        0        0        0 
      50 0.00034403801        0        0        0        0        0 
      60 0.00041365623        0        0        0        0        0 
      70 0.00048208237        0        0        0        0        0 
      80 0.00055122375        0        0        0        0        0 
      90 0.00062060356        0        0        0        0        0 
     100 0.0092129707     3156        0        0        0        0 
     110  0.011072636     3157        0        0        0        0 
     120  0.016091108     3158        0        0        0        0 
     130  0.018477201     3158        0        0        0        0 
     140  0.023477793     3158        0        0        0        0 
     150  0.025398493     3158        0        0        0        0 
     160  0.027075768     3158        0        0        0        0 
     170  0.029419422     3158        0        0        0        0 
     180  0.031110525     3158        0        0        0        0 
     190  0.035495281     3158        0        0        0        0 
     200  0.047222137     6267        0        0        0        0 
     210  0.050610065     6267        0        0        0        0 
     220  0.053887367     6265        0        0        0        0 
     230  0.057132006     6265        0        0        0        0 
     240  0.060388088     6266        0        0        0        0 
     250  0.063811541     6266        0        0        0        0 
     260  0.067081213     6264        0        0        0        0 
     270  0.070362091     6233        0        0        0        0 
     280  0.073603153     6098        0        0        0        0 
     290  0.076751709     5826        0        0        0        0 
     300  0.085217714     8553        0        0        0        0 
     310  0.089679241     8141        0        0        0        0 
     320  0.093842506     7731        0        0        0        0 
     330  0.097784996     7384        0        0        0        0 
     340   0.10156369     7119        0        0        0        0 
     350   0.10542011     6915        0        0        0        0 
     360   0.10906863     6729        0        0        0        0 
     370   0.11265039     6581        0        0        0        0 
     380   0.11612058     6388        0        0        0        0 
     390   0.11947346     6024        0        0        0        0 
     400    0.1281302     8742        0        0        0        0 
     410   0.13270044     8258        0        0        0        0 
     420   0.13695812     7815        0        0        0        0 
     430   0.14104915     7439        0        0        0        0 
     440   0.14494848     7178        0        0        0        0 
     450   0.14892983     6971        0        0        0        0 
     460   0.15258145     6798        0        0        0        0 
     470   0.15622878     6626        0        0        0        0 
     480   0.15970159     6420        0        0        0        0 
     490   0.16300559     6056        0        0        0        0 
     500   0.17162919     8767        0        0        0        0 
     510   0.17622185     8306        0        0        0        0 
     520   0.18047023     7854        0        0        0        0 
     530   0.18449903     7495        0        0        0        0 
     540    0.1883359     7169        0        0        0        0 
     550   0.19234943     6958        0        0        0        0 
     560   0.19608259     6782        0        0        0        0 
     570   0.19960904     6637        0        0        0        0 
     580    0.2030437     6438        0        0        0        0 
     590   0.20635939     6069        0        0        0        0 
     600   0.21502852     8849        0        0        0        0 
     610   0.21965122     8368        0        0        0        0 
     620   0.22393084     7924        0        0        0        0 
     630   0.22797537     7559        0        0        0        0 
     640   0.23183632     7234        0        0        0        0 
     650   0.23575902     6978        0        0        0        0 
     660   0.23938227     6826        0        0        0        0 
     670   0.24292564     6675        0        0        0        0 
     680   0.24638271     6473        0        0        0        0 
     690   0.24971318     6155        0        0        0        0 
     700   0.25844693     8866        0        0        0        0 
     710   0.26307058     8377        0        0        0        0 
     720   0.26737022     7949        0        0        0        0 
     730   0.27141404     7557        0        0        0        0 
     740   0.27527523     7235        0        0        0        0 
     750   0.27920604     7006        0        0        0        0 
     760   0.28282595     6821        0        0        0        0 
     770   0.28637409     6688        0        0        0        0 
     780   0.28982925     6473        0        0        0        0 
     790   0.29316473     6116        0        0        0        0 
     800   0.30183268     8817        0        0        0        0 
     810   0.30655622     8335        0        0        0        0 
     820   0.31082559     7851        0        0        0        0 
     830   0.31482506     7460        0        0        0        0 
     840   0.31866169     7177        0        0        0        0 
     850   0.32256293     6954        0        0        0        0 
     860   0.32617831     6796        0        0        0        0 
     870   0.32971239     6655        0        0        0        0 
     880   0.33316636     6438        0        0        0        0 
     890    0.3364737     6102        0        0        0        0 
     900   0.34511018     8800        0        0        0        0 
     910   0.34971571     8320        0        0        0        0 
     920   0.35397267     7893        0        0        0        0 
     930    0.3579762     7479        0        0        0        0 
     940   0.36180186     7166        0        0        0        0 
     950   0.36569619     6957        0        0        0        0 
     960   0.36933661     6761        0        0        0        0 
     970   0.37283516     6600        0        0        0        0 
     980   0.37624598     6379        0        0        0        0 
     990    0.3795259     6043        0        0        0        0 
    1000    0.3881762     8801        0        0        0        0 
Loop time of 0.388226 on 1 procs for 1000 steps with 8801 particles

MPI task timing breakdown:
Section |  min time  |  avg time  |  max time  |%varavg| %total
---------------------------------------------------------------
Move    | 0.33606    | 0.33606    | 0.33606    |   0.0 | 86.56
Coll    | 0          | 0          | 0          |   0.0 |  0.00
Sort    | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.0018075  | 0.0018075  | 0.0018075  |   0.0 |  0.47
Modify  | 0.040013   | 0.040013   | 0.040013   |   0.0 | 10.31
Output  | 0.0059354  | 0.0059354  | 0.0059354  |   0.0 |  1.53
Other   |            | 0.004408   |            |       |  1.14

Particle moves    = 5939599 (5.94M)
Cells touched     = 6029568 (6.03M)
Particle comms    = 0 (0K)
Boundary collides = 28392 (28.4K)
Boundary exits    = 22819 (22.8K)
SurfColl checks   = 0 (0K)
SurfColl occurs   = 0 (0K)
Surf reactions    = 0 (0K)
Collide attempts  = 0 (0K)
Collide occurs    = 0 (0K)
Reactions         = 0 (0K)
Particles stuck   = 0
Axisymm bad moves = 0

Particle-moves/CPUsec/proc: 1.52993e+07
Particle-moves/step: 5939.6
Cell-touches/particle/step: 1.01515
Particle comm iterations/step: 1
Particle fraction communicated: 0
Particle fraction colliding with boundary: 0.00478012
Particle fraction exiting boundary: 0.00384184
Surface-checks/particle/step: 0
Surface-collisions/particle/step: 0
Surf-reactions/particle/step: 0
Collision-attempts/particle/step: 0
Collisions/particle/step: 0
Reactions/particle/step: 0

Surface reaction tallies:
  id adsorb_test_ps2 style adsorb #-of-reactions 5
    reaction all: 14
    reaction C(b) --> C(g): 14

Particles: 8801 ave 8801 max 8801 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Cells:      8 ave 8 max 8 min
Histogram: 1 0 0 0 0 0 0 0 0 0
GhostCell: 0 ave 0 max 0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
EmptyCell: 0 ave 0 max 0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
